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Abstract 

We consider the mean curvature evolution of rotationally symmetric surfaces. Using 
numerical methods, we detect critical behavior at the threshold of singularity formation 
resembling the one of gravitational collapse. In particular, the mean curvature simulation 
of a one-parameter family of initial data reveals the existence of a critical initial surface 
that develops a degenerate neckpinch. The limiting flow of the Type II singularity is 
accurately modeled by the rotationally symmetric translating soliton. 



1 Introduction 



Geometric evolution equations are diffusive (heat type) equations that describe the de- 
formation of metrics on Riemannian manifolds driven by their curvature in various forms. 
These deformations are naturally divided into two classes that correspond to intrinsic and 
extrinsic curvature flows. The former class refers to deformations driven by the intrinsic 
Ricci curvature tensor on a manifold, whereas the latter refers to deformations of sub- 
manifolds embedded in higher-dimensional spaces evolving by their extrinsic curvature. 
The field of geometric flows is a subject of interest in both physics and mathematics with 
important recent advances allowing for further striking developments. 

The understanding - and in many cases the mathematical classification - of the for- 
mation of possible singularities along different flows has played an important role in the 
progress of the field of geometric flows. The authors of [3J, by numerically simulating 
the Ricci flow of a specific one-parameter family of rotationally symmetric geometries on 
S 3 with a varying amount of S 2 neckpinching, found critical behavior at the threshold 
of singularity formation. According to their results, the Ricci flow for a critical initial 
geometry, i.e., for the transition point between initial geometries on S 3 which under the 
volume normalized flow converge to a sphere and other initial geometries that develop a 
singular neckpinch, evolves into a degenerate neckpinch. In addition, in [J] it was shown 
that the limiting flow of the critical geometry is modeled by the Bryant Ricci soliton. 

In this work we initiate the search for critical behavior of singularity formation of 
extrinsic curvature flows. With the purpose of revealing a critical behavior analogues to 
the one close to the threshold of black hole formation pQ, [2], we perform a numerical 
study of the singularity structure of closed rotationally symmetric surfaces evolving under 
their mean curvature. 

In order to investigate whether critical behavior is present in extrinsic curvature flows, 
we follow a similar numerical approach to the one of [3JH] adjusted to the mean curvature 
flow. The plan of the paper is as follows: section 2 starts by briefly reviewing the mean 
curvature evolution of surfaces of revolution and commenting on the classification of 
their singularity pattern. Then, a one-parameter family of surfaces with a dumbbell- 
like shape is introduced to be used in the current study. Emphasis is placed upon the 
two-dimensional mean curvature translating soliton and its role in Type II singularity 
formation. In section 3, the numerical results are presented. These results provide 
strong evidence that critical behavior is indeed detected in the simulation of the mean 
curvature evolution of the one-parameter family of initial data of our choice. The limiting 
flow of a critical initial surface is shown to be accurately modeled by the appropriate 
mean curvature steady soliton. Furthermore, in the same section, it is shown that the 
critical surface develops a degenerate neckpinch at the final stage of its evolution. Where 
appropriate, frequent references are made to existing mathematical results that formed 
the original basis of the expectations for mean curvature critical behavior. Finally, section 
4 presents both the conclusions and possible future research directions related to critical 
phenomena of gravitational collapse. 
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2 Mean curvature flow 



A smooth n-dimensional hypersurface embedded in Euclidean space, X : M n — > M. n+ , 
evolves by its mean curvature if the normal velocity of any point p G M n , at the time of 
evolution t, coincides with the mean curvature H(p,t) at that point, 

^^ = -H(p,t)n(p,t), (2.1) 

where n is the outer unit normal vector at X(p, In a local coordinate system the 
metric and the second fundamental form induced on M t = X t (M n ) are given respectively 
by 

9af}M = d a X(p,tydpX(p,t) , (2.2) 

and 

hapfat) = -h(p,t)-d 2 Q pX(p,t) , (2.3) 

with 1 < a, (3 < n. The eigenvalues of the second fundamental form h a p with respect 
to the induced metric g a p, usually denoted the principal curvatures of the 

hypersurface M n . The mean curvature is given by the trace of the second fundamental 
form H = g a/3 h a [3 — K\ + k 2 + ■ ■ ■ K n , and is related to the scalar curvature R according 
to 

R(p,t)=H 2 (p,t)-\A(p,t)\\ (2.4) 

where \A\ 2 = g aP g l5 h ai hp S = k\ + k\ + . . . + «£. 

In the following sections, and after first specializing on the particular surfaces of 
interest, the differential equations that govern the evolution of the geometric quantities 
just introduced are used. These are the time evolution equations of g a p, h a p and H (that 
can be found in the introductory chapters of any textbook on mean curvature flow, e.g. 
0), 



dg, 



aii 



-2Hh a(3 , (2.5) 
Ah a(3 - 2Hh ai K 1 p + \A\ 2 h aP , (2.6) 
AH+\A\ 2 H, (2.7) 



dt 

®hgf3 _ Wi , yu h h - | , [-2 / 

dt 

dH 
~dt 

where the dependence of all quantities involved on p and t has been suppressed. 

Under the mean curvature flow (12. ip . a compact and closed hypersurface M of pos- 
itive mean curvature embedded in Euclidean space gets extinct after a maximal finite 
time T* subject to the bounds [HI [7J E] 

o V(M f < < diam ( M °) 2 ; (2 .^ 



A(M { 



1 The signs are chosen such that H = —Hn is the mean curvature vector with the mean curvature of 
a convex hypersurface being positive. 
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given here in terms of the volume V(M ) enclosed by the initial n- dimensional hypersur- 
face, its initial area A(M ) and diameter. In the course of its evolution the hypersurface 
may possibly become singular at an earlier time T < T*. In this case, the solution M t of 
the mean curvature flow (12. ip for the initial hypersurface M exhibits a singularity before 
its collapse to one point (or possibly more), that is classified as Type I or II according 
to the rate at which the maximal curvature A max (t) = m&xA(t) blows up as t — > T. If 

Alt 

the maximal curvature can be shown to obey 

A max (t)<-^==, (2.9) 

for some constant C < oo, then M t exhibits a Type I singularity. Otherwise, it is said 
that a Type II singularity is forming as t —>■ T. Compared to the curve shortening flow, 
note that singularities of the second type, or "slowly forming" singularities, do not occur 
during the evolution of compact curves embedded on the Euclidean plane. The final stage 
of such curves moving under their mean curvature is to become convex - if they were not 
initially - and shrink to a point as circles [9]. Similarly, compact convex hypersurfaces 
asymptotically converge to the self-similar round sphere solution and collapse to a point 
of K n+1 with the formation of a Type I singularity |10j . 

The prototype example of a flow that is commonly invoked to argue in support of 
the possibility that an embedded hypersurface may become singular before it shrinks to 
a point is provided by an appropriate dumbbell in M. 3 , see for instance jTTj. Consider a 
dumbbell consisting of two large approximately round spheres connected by a sufficiently 
thin cylindrical neck as initial configuration. Due to its large curvature, the neck pitches 
and a singularity is formed before the two spheres have the chance to shrink considerably. 
However, if the neck is not "very thin", the dumbbell surface becomes convex in the 
course of its flow and, in accordance to the standard theorem of (TUj, collapses to a single 
point. An interesting intermediate scenario is also possible. There exist dumbbell- shaped 
surfaces with just the right neck that shrink to a point without ever turning convex nor 
loosing their neck under their mean curvature evolution. This specific class of mean 
curvature flow solutions that exhibit precisely this degenerate neckpinching behavior is 
the primary interest of the current work. The maximal curvature for these solutions does 
not satisfy (12.91) for any finite constant C, thus the forming singularity at their collapse is 
of Type II. The existence of such mean curvature flow solutions was originally suggested 
by R. Hamilton [12]. He indicated that solutions with a maximal curvature blowing up 
faster than the rate ( 12. 91) can be produced by studying rotationally symmetric surfaces 
that exhibit non-generic neckpinching. These solutions were investigated using level-set 
methods initially in [7J and later in [12] where a proof of the Type II character of the 
forming singularity was presented. A detailed study of the possible blow-ups of these 
solutions can be found in [13] along with a piecewise construction of the corresponding 
solutions. 

It should be noted here that the above discussion refers to left-right symmetric 
dumbbell-shaped surfaces, in accordance with the choice of surfaces to be investigated in 
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the sequel. Reflection non-symmetric surfaces with an intermediate, or "critical" in the 
nomenclature of the present work, initial data may also exhibit similar singular behavior 
[T2| [13]. This discussion can be easily generalized to hypersurfaces of dimension n > 3 
and/or more than one number of initial necks [T3j. As a final point, a reference is 
made to p3] and [15] for the curve shortening evolution and Type II singularity formation 
of planar or space not embedded curves. 



2.1 Surfaces of revolution 

In this work rotationally symmetric compact surfaces of a dumbbell shape are employed as 
initial configurations to study the formation of Type II mean curvature flow singularities. 
The expression that a generic dumbbell- shaped surface exhibits a pinch at some time T < 
T* (see equation (I2.8P ) is used if, as t — > T, its neck pinches off leaving behind two "tear" 
shaped components on either side of the neck. This work will not be concerned with the 
subsequent evolution of these two separate components. For a rigorous presentation of 
the generalized solution beyond the first singular time T using level-set flow methods, 
see e.g. [H [121 EE E] an d for an update on surgery techniques see [18j. Next, the 
essentials for describing the mean curvature evolution of a compact surface of revolution 
are introduced and an adequate as well as simple set-up is provided to accomplish the 
task of numerically revealing critical behavior in mean curvature flow. 

A surface of revolution is produced by revolving a planar profile curve (generator) 
about a line in K 3 (the axis of revolution). For convenience, the profile curve is chosen 
to lie in the xy-plane of M. 3 and the axis of revolution to coincide with the x-axis. Using 
the principal patch of the standard parametrization of a compact surface of revolution 
generated by a curve (x(u),y(u)), 

X = x{u) , Y = y{u) cos0 , Z = y{u) sin0 , (2-10) 

with <p G (0, 27r) and u ranging in some finite interval, the induced area element has 
components 

g uu = x\uf + y'(u) 2 , g H = y 2 (u), ^ = 0, (2.11) 
while the unit surface normal is 

h = — , (y'(u), —x'(u) cos0, —x'(u) sin</>) . (2.12) 

The mean curvature H (u) of a surface of revolution is given by the sum of the principal 
curvature of its meridians k u and the principal curvature of its parallels 

x'(u)y"(u)-y'(u)x"(u) 
(x'(u) z + y (u) z ) 1 
x'(u) 

H = — m ■ 2 - 14 
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The time evolution of a rotationally symmetric surface embedded in R 3 under its mean 
curvature 

^^ = -H(u,t)n(u,t), (2.15) 

can effectively be encoded into the time evolution of the revolving planar curve that 
generates it. This way, the mean curvature flow in R 3 is reduced to a planar deformation 
of a certain kind for the profile curve described as a graph form in the following subsection. 

Under mean curvature flow, an initially rotationally symmetric surface Mq remains 
rotationally symmetric in the course of its regular evolution M t . The evolution of the 
induced metric is dictated by (12.5ft in terms of the mean curvature and the components 
of the second fundamental form. More precisely, if the non-zero components of the two- 
dimensional metric are denoted as 

g uu (u,t) = S ,2 (u,t) + R ,2 (u,t), (2.16) 
g H {u,t) ee R 2 {u,t) , (2.17) 

then it can easily be verified that the flow equations ( 12.51) for the metric of a surface of 
revolution can be described by the planar evolution 

OR S' S' 2 

— — = n (S R — R S ) t — o sr , (2.18) 

dt (S' 2 + R' 2 ) R(S' 2 + R' 2 ) v ; 

= o (R'S" - S'R") + . 9 ^ , (2.19) 

dt (S> 2 + R' 2 ) R(S' 2 + R' 2 ) ' v ; 

for the two Cartesian coordinates R(u,t) = y(u,t) and S(u,t) = x(u,t) of the profile 
curve that generates it. 

In contrast to the case of the formation of Type I singularities, and to the best of our 
knowledge, there exist no mean curvature flow solutions know in closed form that exhibit 
a Type II singularity. This is also the case for the curve shortening flow of even planar 
curves for which, apart from the homothetically contracting circle and the self-similar 
solutions with intersections of [19] , the set of compact known solutions is exhausted by 
the oval-shaped solutions that appear in [TT] (paper-clip model [2Ql EI]). The latter 
class of solutions can be understood as the mean curvature flow analogue of the Ricci 
flow Rosenau solution [22J (or sausage model [23]). At its final stage, this compact 
embedded curve flowing under its mean curvature collapses to a single point while a 
Type I singularity is formed. Due to the lack of existence of ancient mean curvature 
flow solutions which can be followed until the occurrence of a Type II singularity, it is 
necessary to turn to numerical methods to pursue the investigation of the formation of a 
Type II singularity during the evolution of a compact surface of revolution. The initial 
data is provided by the surfaces produced by the revolution of a particular family of 
closed planar curves. The study can nevertheless be extended by incorporating different 
families of surfaces with similar asymptotic behavior. 
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The Cassini ovals are planar quartic curves that constitute a generalization of the 
lemniscate of Bernoulli, and in that respect of the ellipse. A Cassini curve is described 
as the locus of points for which the product of the distances to each of two focal points 
separated by a distance 2a equals a constant b 2 . The corresponding Cartesian equation 
is 

(x 2 + i/ 2 + a 2 ) 2 -4aV-6 4 = 0. 

The shape of the Cassini ovals depends on the ratio A = | which will be denoted as shape 
parameter from now on. Despite their name, these curves do not have an oval shape for 
generic values of the shape parameter. For A < 1 the curve is a single loop with an oval 
or peanut shape, for A = 1 the lemniscate curve is produced, while for A > 1 they are 
comprised of two disjointed loops. 




A = A < 1 A = l A > 1 



Figure 1: Cassini curves for different values of the shape parameter A. 

The smooth surfaces of revolution produced by revolving the family of the Cassini 
curves of a single loop with A < 1 around their horizontal axis (Figured]), provide us with 
the initial data M . Based on equations (12.181) and (12.191) . the mean curvature evolution 
M t is followed without concerning if it is possible to be followed backwards in time to some 
ancient solution. It should be noted that the resulting rotationally symmetric surfaces 
are convex for A < 4g and, as a consequence, the surfaces M t evolving under their mean 
curvature unavoidably shrink to a point by acquiring a spherical shape, within this range 
of A. On the other hand, the surfaces generated by Cassini curves with ^= < A < 1, 
that constitute our variant of a dumbbell surface, may possibly pinch depending on the 
precise value of the shape parameter A. The detailed numerical evolution for our set-up is 
presented in section 3, where it is shown that, for a critical value of the shape parameter, 
the flow terminates with the occurrence of a degenerate neckpinch. 



2.2 Effective one-dimensional flow and the translating soliton 

A solution of the mean curvature flow (12. ip that exists for all time and moves by trans- 
lation in M n+1 is called a translating soliton. More generally, the interest in solitonic 
solutions of various geometric flows is, among other reasons, due to the fact that such 
solutions model the asymptotic behavior of developing singularities. Regarding the for- 
mation of Type II singularities, recall the Bryant soliton as a limiting flow of a three- 
dimensional Ricci flow solution exhibiting such a singularity; or the emergence of the 
grim-reaper solution during the analogue scenario in the curve shortening flow [HI 115 J 2 ! 

2 For a treatment along these lines in relation to the two-dimensional cigar Ricci soliton see [24] . 
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The study of hypersurfaces moving under their mean curvature near a singularity us- 
ing rescaling techniques [23] involves solitonic solutions. The relevant two-dimensional 
translational solution for the mean curvature flow has been studied among other people 
by the authors of 0, EH [E0 [25] , with early numerical evidence for its existence that can 
be traced in [26]. In [27] it was proven that there exists a convex, rotationally symmetric 
solution translating under its mean curvature with any prescribed constant speed. Un- 
fortunately, an explicit solution describing this soliton in closed form is not available in 
the literature. 

Referring back to the set-up of section 2.1, consider a planar evolving curve T t (x,y) 
that can be represented either as a horizontal graph y = y(x,t), or as a vertical graph 
x = x(y, t). As long as the evolution M t of the surface of revolution generated by rotating 
r (x, y) around the x-axis is regular, the mean curvature flow evolution (12.151) in R 3 is 
dimensionally reduced to the planar deformation for the graph of the underlying curve 
T t (x,y) described either by the horizontal modified, in reference to the curve shortening 
flow, graph equation 

dy(x,t) y"(x,t) 1 



dt " l + y' 2 (x,t) y(x,t)' 
or equivalently by the corresponding vertical graph equation 



(2.20) 



dx(y,t) x"(y,t) x'(y,t) 

dt l+x' 2 (y,t) + y ' 1 ' 

where prime denotes differentiation with respect to the corresponding spatial variable. 
The last term of equations (12.201) and (12.211) captures the S 1 extrinsic curvature of the 
surface due to the revolution around the x-axis, and differentiates the resulting one- 
dimensional flow from the ordinary curve shortening flow for the graph of its generator. 

Within this effective description of the mean curvature evolution of a rotationally 
symmetric surface by means of the flow equation of its profile curve (I2.20p or (I2.2ip . 
a suitable presentation of the corresponding two-dimensional translating solitonic solu- 
tion follows by considering the possible reparametrizations along this flow. By properly 
introducing reparametrizations generated by a vector field £ to take into account the 
translational motion of the soliton (see for instance [20]), the mean curvature flow in R 3 
( I2.15P acquires the form 

-^ = -Hh + l (2.22) 

This way, a vector field £ = df(x,y) captures the rigid motion of a gradient soliton 
that translates with a constant velocity c > along the positive x-axis for f(x) = ex. 
Denoting it as a horizontal graph y = <p(x), one easily verifies that it satisfies the ordinary 
differential equation 

V {x) V "{x) = (1 + V '\x)) (1 - c<p(x)<f/(x)) , (2.23) 

with respect to the abscissa of the planar curve x, or its equivalent one following from 
( I2.21j) when is preferable to consider a vertical graph. A numerical plot of the solution 
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of (12.231) is presented in Figure along with a 3D graphical illustration of the two- 
dimensional mean curvature flow translating soliton. Alternatively, by canceling the 
"freezing" effect of the reparametrization term in (I2.22p . the generator curve of the soliton 
solution may be viewed as moving linearly in time along the x-direction with constant 
speed c, 

y(x,t)=ip(x + ct). (2.24) 

Likewise, the profile curve that generates a translating solution along the ?/-axis, arises 
by reversing the roles of the coordinates x and y and considering vertical graphs instead. 
Furthermore, the sign of the velocity c selects the orientation of the translational solution 
in both cases, distinguishing a soliton from an anti-soliton solution. In general, the 
rotationally symmetric gradient soliton that is generated in this way, with any line of R 3 
as axis of symmetry along which it is translating with a non-zero finite velocity c, has a 
"tip" at its intercept with that line and asymptotically opens to a hyperboloid. 





(a) 



Figure 2: The planar curve (a) that generates the "bowl" soliton (b) translating with 
unit speed to the right. 



A convenient way to describe the mean curvature of the soliton under investigation 
is to introduce the slope of the curve that produces it via rotation. In general, the mean 
curvature if of a smooth rotationally symmetric surface generated by revolving a profile 
planar curve y(x) around the x-axis is related to the mean curvature of its generator 



curve H c as 



H = H C + 



cos/? 

y 



(2.25) 



where j3 is the angle formed by the tangent at each point of the curve and the x-axis, 
i.e., 

P = arctany'(x) , (2.26) 

when viewed as a horizontal graph. In terms of this angle, the mean curvature of the 
soliton translating along the positive x-axis with velocity c acquires the simple form 

H sol = c sin/3 with 0</?<vr, (2.27) 



where the defining equation (12. 23[) was used. It is instructive to realize that the result 
( 12.27!) coincides with the mean curvature of the translating soliton of the curve shortening 
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flow on the two-dimensional Euclidean plane, the grim-reaper solution [20J. Put differ- 
ently, the planar curve that satisfies (I2.23P is such that the rotation around the x-axis 
produces just the curvature required so that the generated non-compact soliton surface 
has mean curvature, in terms of its slope, matching the one of its one-dimensional cousin. 
The expression (I2.27P for the mean curvature of the two-dimensional soliton, given here 
in closed form in terms of the slope with respect to its axis of symmetry, corresponds to a 
static solution of equation (12.71) when reparametrizations are taken into account. Equa- 
tion (12.71) may in principle be used to study the stability properties of the rotationally 
invariant translating soliton solution [28], along the lines of [20J. For completeness, in 
Figure [3J the mean curvature of the soliton is presented as a function of the horizontal 
distance from its tip. The inspection of this plot reveals that for the soliton (12.231) . the 



sol io 
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Figure 3: Mean curvature of the unit speed translating soliton in terms of the horizontal 
distance from its tip. 



mean curvature is everywhere positive, acquires its maximum value at the tip {(5 = 7r/2) 
and monotonically dies off to zero (/3 — > 0). This is in full accordance with a general 
theorem stating that any eternal strictly convex solution of the mean curvature flow with 
the mean curvature attaining a maximum value at only one point, must necessarily be 
a translating soliton [29]. This is an important ingredient that goes into the study of 
Type II singularities by means of rescaling techniques with the soliton solutions coming 
into play as limiting flows, which will be discussed below. In passing, we note that mod- 
ulo direction, orientation and translational speed, the soliton depicted in Figure [2] is the 
unique two-dimensional translating soliton solution that can be generated by rotating 
a horizontal with respect to the axis of revolution graph. It must be stressed though, 
that other two-dimensional mean curvature flow translating solitons, not falling into this 
category, do existfE 



3 We thank G. Huisken for pointing out their existence and also clarifying their possible relation to 
the formation of Type II singularities in the course of the mean curvature flow evolution of immersed 
(hyper)surfaces (see also [30]). 
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2.3 Rescaling the singularity 



The asymptotic description of a mean curvature flow solution near a singularity of both 
types, through the examination of sequences of rescaled solutions, involves the notion of 
the so-called blow-ups. Consider the evolution M t of a compact, smooth n-dimensional 
hypersurface of positive mean curvature for t G [0, T), where T is the first singular time. 
Let po be the point of Mt at which the curvature becomes unbounded as t —>■ T (the 
blow-up point), and the forming singularity be of Type II. If for every integer k > 1, 
t k G [0, T — |] and p k G M t , {p k , t k } is a sequence such that 

H 2 (p k ,t k )(T-j-t k ) = max H 2 (p, t)(T - \ - 1) (2.28) 

k t<T-i/k,peM t k 

then the sequence of the rescaled flows of X : M t — > R n+1 , 

X {k) (p,r) =t k (x{p,t k + tfT)-X{p k ,t k )) for re[-4t k ,e 2 k (T-t k -^)}, (2.29) 
where 

6 k = H(p k7 t k ), (2.30) 

converges smoothly to a mean curvature eternal flow. The limiting hypersurface is weakly 
convex and has uniformly bounded mean curvature [8]. Section 3 provides numerical 
evidence that the type II blow-up fl 2 . 2 9 [) of the mean curvature flow of the dumbbell- 
shaped surface of revolution generated by the Cassini curve with shape parameter equal 
to the critical value A c , with either of its two poles as blow-up point, converges to the 
rotationally symmetric translating soliton (12.23)) . 

In the same section, there is also a focus on the limiting flow of the central region of 
the neck of the evolving surfaces where the pinching is observed for generic supercritical 
values of A. In [25] (but see also [31]) it was shown that the limiting flow of the parabolic, 
or type I, blow-up of a rotationally symmetric shrinking neck converges to the cylinder 
of radius y/2. In terms of the profile curve that generates such a surface by rotation, this 
amounts to 

; 1 y ( ; X ) -> V2 as t->T. (2.31) 

The continuous rescaling (12.311) is relevant for the asymptotic description of the generic 
neckpinch that occurs when the regular evolution of our supercritical surfaces comes to an 
end with the formation of a Type I singularity. For the case of a degenerate neckpinching 
a more refined rescaling is needed. This is provided in [13], and in terms of the profile 
curve it has the form 

y(x) = V2 + KB.m m (x){T -t)^- 1 + 0{{T -t)^ 1 ) as t->T, (2.32) 

with 

y=(T-ty 1/2 y and x = (T-t) 1/2 x, (2.33) 

where Hm m (i) is m th Hermite polynomial normalized such that the x m term has unit 
coefficient. For a left-right symmetric surfaces of revolution, the integer m appearing in 
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(I2.32p is even and greater than two. Also, in that case, the numerical coefficient K is 
determined based on the distance between the Type II singular point and the poles of 
the surface. 

3 Numerical investigation 

In this section the results of the numerical simulations of the mean curvature flow of 
the chosen one-parameter family of rotationally symmetric dumbbell surfaces in R 3 are 
presented. The initial configurations are provided by the compact surfaces of revolution 
generated by revolving the one-component Cassini curves around their horizontal axis 
of symmetry, as described in section 2. These are parametrized in terms of the shape 
parameter, with < A < 1. For small values of this parameter ("loose corsetting" relative 
to a round sphere in the terminology of 0, H]), even if the initial data is not convex, the 
surfaces evolving under their mean curvature become convex and asymptotically approach 
the geometry of the round two-sphere. In contrast, the dumbbell surfaces with large 
values of the shape parameter A ("tight corsetting") do not loose their neck until their 
evolution becomes singular with the development of a generic neckpinch at the center 
of their left-right symmetric extend. The numerical simulation of the mean curvature 




Figure 4: A generic dumbbell-shaped supercritical initial surface. 

flow of such initial data with two kinds of possible asymptotic behavior depending on 
the shape parameter A, reveals the existence of a critical transition value A c ~ 0.9076 for 
which a degenerate neckpinch occurs. Unlike both subcritical (A < A c ) and supercritical 
(A > A c ) initial surfaces, the smooth evolution of the dumbbell with A = A c terminates 
with the formation of a Type II singularity. Focusing on the flow of this critical initial 
data, it is possible to verify that the critical surface in the course of its evolution becomes 
increasingly cylindrical in the neighborhood of the initial pinching and that the mean 
curvature blows up at its two poles. This is in full analogy to the critical flow of [3j 
H] for the Ricci flow. It is valuable to stress that, at its very final stage, the critical 
dumbbell surface collapses to a point but without turning convex until that time. The 
close examination of the limiting flow, with either of the two poles as a blow-up point, 
provides strong numerical evidence that it is modeled with high accuracy by the two- 
dimensional translational mean curvature soliton (I2.23p . 

The mean curvature evolution of this one-parameter family of surfaces is studied by 
numerical means. One of the parameters of the Cassini closed curves that generate them 
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is held fixed, b = 1, and the other one takes values in the interval [0, 1). This way, the 
induced metric of the initial data is parametrized in terms of the scale parameter A, 

ds 2 = T 2 {9 ■ A) d9 2 + R 2 {6 ■ A) d<f) 2 , (3.34) 

where T 2 (9;\) = >^ 2 (#;A) + R' 2 ( y 9;\) with < A < 1 and < 9 < it. The metric 
components are numerically evolved based on the mean curvature evolution equations 
(I2.18P and (12.191) for the two planar coordinates S(9, t) and R(9, t) of the profile curve as 
functions of the angle 9 and the evolution time t. Angular derivatives are approximated 
by centered finite differences with an angular step size A9 = n/N, where N is the 
number of angular grid points 9i employed. As in 0, H] the angular grid is augmented 
by considering two external fictitious points 9q and 9n+i on both of its ends. Periodic 
boundary conditions are implemented by imposing the conditions S (t = 0) = Si(0), 
Sat + i(0) = Sn(0) and -Ro(O) = — -Ri(O), Rn+i(0) = — -Rjv(O) at the two poles. In order 
to meet the need of increased accuracy close to the poles of the evolving surface {9 = 
and 9 = n), especially when approaching the forming singularity in time, variable 
temporal step size is used when the differentiation with respect to time in (12.181) and 
(I2.19p is approximated by finite differences. Computer code was executed in Mathematica 
with up to 10,000 angular grid points. This was reduced by a factor of two when the 
shape parameter of the initial surface was adjusted close to its critical value in order 
to avoid numerical instabilities due to the stiffness of the simulated system of finite 
difference equations. In general, the results of repeated runs of the computer code clearly 
differentiate the asymptotic behavior of subcritical and supercritical initial geometries in 
a well distinguishable way. 

3.1 Critical flow and comparison with the soliton 

As in pE], two complementary comparison methods are presented in order to confirm the 
claim that the flow of the critical initial surface is asymptotically modeled by the mean 
curvature translating soliton. 

As a first examination approach, the geometry of an appropriate mean curvature 
soliton solution is compared with the geometries near the left pole for a sequence of 
supercritical initial geometries, with shape parameters approaching the critical value 
A c , at the time of maximum curvature. The soliton's speed parameter c (cf. section 
2.2) is determined by requiring in each case the mean curvature at its tip to match the 
maximum value of the mean curvature at the poles of the evolving dumbbell surfaces. 
Note that during the mean curvature evolution of these surfaces, the maximum value of 
the curvature at the poles is attained when a generic neckpinch occurs at their center 
rendering any evolution thereafter singular. A comparison is presented in Figure [5j where 
the profile curve that generates the corresponding soliton (solid line) is plotted against the 
profile curve that, by revolution, produces the final phase of three different supercritical 
initial geometries while their neck is pinching (dashed lines) versus the horizontal distance 
from their left pole up to their center. Due to the unequal spatial extend of the dumbbells 
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at this singular time, three separate graphs are depicted. Note that the translational 
speed parameter of the soliton increases from left to right. By inspection of this set of 




Figure 5: The profile curve for three supercritical surfaces at pinching (dashed lines) 
compared to the soliton curve (solid lines). Only the part from the left pole up to the 
center is depicted and Ai = 0.94, A 2 = 0.93, A 3 = 0.92. 

graphs, it is clear that as the shape parameter A of the initial data decreases towards the 
critical value, the surface near the pole of the final stage of the singular dumbbell tends to 
approach the surface of the corresponding mean curvature soliton solution. Recall that, 
in this comparison test, no magnification of the dumbbell's surface has been performed 
during its evolution. The accuracy of the polar modeling of the dumbbell by the soliton 
increases for initial data corresponding to slightly supercritical shape parameter values. 
Evidently, owing to the reflection symmetry of the family of the initial configurations, a 
similar behavior is observed when the flow close to the right pole is examined. 

In order to obtain an accurate description of the critical flow corresponding to an 
initial data with a critical shape parameter value A c , a second comparison analysis is 
performed. This time, a choice for a critical initial surface is made with shape parameter 
as close to its critical value as permitted by the numerical determination of A c . The 
flow of the dumbbell with A ~ A c is simulated and it is observed that the maximum 
curvature at the poles increases along the flow of this close to critical initial surface. It 
eventually diverges when I is fine-tuned, as attainable, to be critical. Next, a sequence 
of times approaching the time of singularity formation is considered. For each chosen 
instant of evolution time, the blow-up geometry as specified in section 2.3 is determined, 
with one of the two poles as blow-up point. This way the rescaled dumbbell surface has 
maximal curvature at the left, for example, pole equal to one. This is repeated for the 
remaining subsequent time instances of the critical flow. The results of this procedure are 
graphically demonstrated in the following figures. In Figure EJ the blow-up of the profile 
curve for six different times of the evolution of the "critical" initial surface are plotted 
in comparison with the profile curve that generates the mean curvature steady soliton 
(solid line) corresponding to a translational speed of unit value to the right, in terms of 
the horizontal distance from the tip of the soliton. It is observed that the rescaled curves 
near the pole approach the profile curve of the unit speed soliton solution with increasing 
agreement along the flow of the critical surface. The presentation of the results of this 
simulation is supplemented by the graphical comparison of the mean curvature near the 
poles for a sequence of rescalings of the critical flow to the mean curvature close to the tip 
of the steady soliton of Figure [3 The simulated mean curvature for subsequent instances 
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Figure 6: A sequence of subsequent rescaled profile (dashed) curves along the critical flow 
with time increasing upwards compared to the unit speed soliton profile (solid) curve. 

of time is plotted as a function of the radial distance and excellent agreement with the 
mean curvature of the soliton is detected at later times. 
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Figure 7: Mean curvature near the pole of the rescaled critical surface for six subsequent 
times ti compared to the mean curvature of the soliton (single solid line). 

The results of both of the above described comparison tests provide convincing evi- 
dence that the limiting flow of the critical initial geometry lying at the "transition point" 
of our one-parameter family of initial surfaces, with different mean curvature asymptotic 
behavior on both sides, is accurately modeled by the translating soliton of section 2. 

3.2 Degenerate neckp inching 

The accurate modeling of the polar flow by the surface of the appropriate gradient soliton 
strongly supports the identification of the singularity along the critical flow as a Type 
II. A supplementary verification may also be produced through the detailed study of the 
limiting flow, not close to the poles, but close to the center of the evolving critical surface 
instead. The occurrence of a mean curvature flow degenerate neckpinch can indeed be 
confirmed through its asymptotic description. 
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Certainly, the decisive factor of the classification of a mean curvature singularity is 
the blow-up rate of the maximal mean curvature. Based on the results of [13], the mean 
curvature of the singular point of a rotationally symmetric surface has a limiting form 
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^ + 0(1)) (T-t) 
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1/m-l 



as t 



(3.35) 



where for the case of a left-right symmetric surface m is a positive even integer and £ 
a constant determined by the distance of the poles of the surface from its center at any 
instant of time close to the singular time T. In the case of a Type II singularity, the 
divergence of the mean curvature at the pole of a surface of revolution with one neck is 
expected to be reproduced by the expression (13.351) for m = 4. On the other hand, a Type 
I singularity corresponds to m = 2 in agreement with the rate (12. 9p . Referring to the 
numerical simulations, recall that the evolution of all supercritical initial data becomes 
singular with the occurrence of a generic neckpinch at their center (x — 0) . According to 
[T3] (section 9), the asymptotic shape of the curve that produces this generic neckpinch 
by rotation, while a Type I singularity is formed, is described by 



y{x,T)=K 



\x\ 



log 1 



(3.36) 



where again the constant K is determined by the distance of the poles from the center 
and T is the singular time. In Figure El a graphical comparison of the limiting shape of 
the central region of the profile curve for a typical supercritical surface and the generic 
pinching (13.361) is presented. Adequate agreement is observed as x — > 0. Moreover, the 
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Figure 8: The cusp singularity of a generic neckpinch for A = 0.96; one every 100 grid 
points is depicted out of 10,000 used in total for the simulation. On the right the mean 
curvature of the center point (solid line) is compared to the asymptotic expression (I3.35p . 

Type I characterization of this singularity is also confirmed by the second graph of the 
same figure, where the mean curvature of the central point is compared to the limiting 
blow-up rate (13. 35p with m = 2. In both graphs, the comparison is extended beyond the 
expected agreement. 

A direct numerical verification of the formation of a Type II singularity at the final 
stage of the critical surface based on the expected blow-up rate (13.351) is essentially 
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impossible for two reasons. First of all, for the one-parameter family of initial geometries 
selected here, a degenerate neckpinching occurs for a single value of the shape parameter 
A. As a consequence, the numerical distinction between different blow-up rates (13.351) 
with m = 2 and m = 4 requires an input of very high precision for the critical value 
A c . Second, numerical implications due to the left-right symmetry of the surfaces used in 
the present set-up also prevent such a resolution. Specifically, the coefficient £ in (13.351) 
acquires a vanishing limiting value, since the critical surface collapses to a point as the 
singularity of the second type is formed. 

However, numerical evidence in favor of the occurrence of a degenerate neckpinch 
along the critical flow may be produced by examining the asymptotic shape of the central 
region of the critical surface. This is achieved by considering a sequence of subcritical 
initial surfaces with shape parameter close to the critical value and focusing on their 
central region. During their mean curvature evolution these surfaces remain non-convex 
until a time just before their collapse. The fine examination of the asymptotic central 
shape of these slightly subcritical surfaces just before they turn convex and shrink to a 
point soon after, reveals a match to the degenerate neckpinching of equation (12.321) . In 




Figure 9: A simulation of the degenerate neckpinch for A = 0.9076 with one every 50 
grid points depicted out of a total 5,000 in comparison with the blow-up of the curve 

Figure [HI the central part of the curve that generates such a surface of revolution close to 
its extinction time T is compared to the blow-up of the degenerate neckpinch (solid line) 
described by equation (I2.32H . This plot depicts the generator curve at t = 0.99T (i.e., a 
vertical magnification factor of the order « 50) with initial shape parameter A = 0.9076. 
In this case, improved agreement is detected as the critical shape parameter is approached 
from below, to the extend permitted by the numerically determined value of A c , and at 
later time. This supports the claim that the critical initial surface evolves at its final 
stage into a degenerate neckpinch. Actually, this reasoning can be inverted and used to 
improve the numerical value of A c by demanding accurate modeling of the limiting form 
of the central region by equation (12.321) for the critical surface at the extinction time. A 
precise determination of the critical value of the shape parameter based on non-numerical 
alternative methods would in general be very helpful. 

This presentation of the numerical results concludes with a graphical representation 
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of the dependence of the mean curvature maximum value at the poles of the family 
of initial surfaces on the shape parameter A. For all supercritical surfaces, the mean 
curvature at the poles attains its maximum value at the occurrence of a generic pinching 
at their center. These values are depicted in Figure [TU] for the supercritical range of the 
shape parameter. It is observed that the mean curvature maximum value increases as the 
shape parameter decreases. The depicted divergence as the critical value is approached 
from above is attributed to the Type II singularity of the degenerate neckpinch. The 
numerically extracted values are reproduced by an expression of the form 

H polc = A (A-A c )" n , (3.37) 

where Ao is a constant that depends on the choice of the one-parameter family of initial 
surfaces. For the surfaces of revolution generated by the Cassini curves, the value of the 
mean curvature at the poles of the surface of revolution generated by the lemniscate curve 
(A = 1) may be used as an input. The simulation of the mean curvature flow of this family 
of initial surfaces parameterized by A suggest a value n ~ 1.26 for the "critical exponent" 
of (13.371) . It must be stressed though, that this numerical value highly correlates with 
the value used for A r . 
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Figure 10: Maximal mean curvature of the poles as a function of the shape parameter. 



4 Conclusions and outlook 

In this work we have performed a numerical study of the mean curvature evolution of 
two-dimensional surfaces of revolution of a dumbbell shap^l We considered a specific 
one-parameter family of initial reflection symmetric surfaces generated by revolving the 
Cassini curves around their axis of symmetry. By simulating the mean curvature flow of 

For an early computer study of a collapsing dumbbell under mean curvature flow, that motivated 
investigations of developing singularities before a generic embedded surface shrinks to a point [9], see 
[32]. 
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this family of initial data, parametrized in terms of the shape parameter A of the Cassini 
curves, we succeeded in revealing a behavior in the emerging singularity pattern that 
resembles the critical behavior of gravitational collapse. In particular, all initial surfaces 
with a shape parameter less than a critical value l c collapse to a point in finite time 
by acquiring spherical shape. On the other hand, the evolution of the initial dumbbell 
surfaces with supercritical shape parameter becomes singular with the formation of a 
generic neckpinch at their center. The fine examination of the flow of the critical surface 
allows for the detection of the occurrence of a Type II singularity. Compared to the 
approach of [21 H] on Ricci flow critical behavior, we note that in the present work a 
normalized variant of the mean curvature flow was not employed. The area preserving 
flow of [25], or the modified flow [33J that preserves the enclosed volume, could have 
been incorporated into our study without altering its results. The use of the ordinary 
(unnormalized) mean curvature flow permitted a direct comparison with the results of 
[T3] on Type II mean curvature singularities and allowed for a numerical verification 
of the formation of a degenerate neckpinch at the final stage of the flow of the critical 
surface. 

The numerical results of the mean curvature simulation of the rotationally symmetric 
family of our initial surfaces provide strong evidence that the polar critical flow is accu- 
rately modeled by the rotationally invariant translating soliton solution (12. 23ft . We are 
able to extend this study by incorporating the simulation of different families of initial 
compact rotationally symmetric surfaces, including the one used in [3] and [IJ. In all cases 
examined, mean curvature flow critical behavior is in general observed together with the 
occurrence of a degenerate neckpinch. Critical behavior is also detected by direct numer- 
ical simulations of initial configurations which are not left-right symmetric. In this case, 
the part of the surface on one side of the degenerately collapsing neck shrinks to zero 
size while a Type II singularity is formed. In addition, a similar behavior is observed for 
initial data with more than one neck without any clear evidence of coalescence. Further 
studies are nevertheless required before any universality claims a la Choptuik can be 
made about critical behavior at the threshold of mean curvature singularity formation. 

A relevant future direction of investigation would be to consider the evolution of non- 
rotationally symmetric surfaces within the framework of the current study. This is a very 
challenging and difficult task. As an alternative direction, we suggest the investigation of 
whether critical behavior is present in the mean curvature evolution of torus-like initial 
configurations, as opposed to the dumbbell-like examined in this work ([31] may turn 
out to be helpful in this respect). If critical behavior is detected for this alternative 
class of initial mean curvature data, the comparison of the corresponding critical limiting 
flow with the one of the degenerately collapsing critical dumbbell of the present work 
may reveal the mean curvature analogue of the scale coordinate of gravitational critical 
phenomena [U |2] . 

Mean curvature flow is the basic example of an extrinsic curvature flow. Critical 
behavior and similarities with Choptuik scaling can be sought for other extrinsic flows. 
Taking into account the important role of the inverse mean curvature flow in the proof of 
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the Riemannian analogue of the Penrose inequality of general relativity [31], we believe 
the investigation of critical behavior of this flow deserves an independent study. Such 
a study can motivate the introduction of the appropriate functionals with power-law 
scaling even for the mean curvature flow and may also provide a direct link between 
critical behavior at the threshold of singularity formation of curvature flows and critical 
behavior at the threshold of black hole creation. 

An essential ingredient of the presentation of our study was the translating mean 
curvature gradient soliton that deserves further consideration. Anticipating the role of 
the grim-reaper (hair-pin) solution, as the curve shortening flow translating soliton, in 
the construction of integrable quantum field theories of boundary interactions |21j, and 
also, the role of translating solitons in the level-set approach of the mean curvature flow 
[TT] we indicate that work is in progress towards finding the proper placement of this 
solitonic surface within the context of boundary conformal field theory. Finally, also 
interesting is the investigation of generalizations of this solitonic solution for non-flat 
ambient spaces [20j and/or in the presence of anti-symmetric tensor fields, motivated by 
string theory [351 [36] , that may influence the structure of singularity formation. 
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